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Abstract 

In this paper we investigate the possible direct, non-gravitational interaction between holographic dark energy 
(HDE) and dark matter. Firstly, we start with two simple models with the interaction terms Q oc p^,,, and Q oc p^^., and 
then we move on to the general form Q oc pf„ffj^. The cosmological constraints of the models are obtained from the 
joint analysis of the present Union2.1+BAO+CMB+i/o data. We find that the data slightly favor an energy flow from 
dark matter to dark energy, although the original HDE model still lies in the 95.4% confidence level (CL) region. For 
all models we find c < 1 at the 95.4% CL. We show that compared with the cosmic expansion, the eff'ect of interaction 
on the evolution of pj„, and pde is smaller, and the relative increment (decrement) amount of the energy in the dark 
matter component is constrained to be less than 9% (15%) at the 95.4% CL. By introducing the interaction, we find 
that even when c < 1 the big rip still can be avoided due to the existence of a de Sitter solution at z — » - L We show 
that this solution can not be accomplished in the two simple models, while for the general model such a solution can 
be achieved with a large and the big rip may be avoided at the 95.4% CL. 
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I. INTRODUCTION 



Cosmological observations such as the type la supemovae (SNIa) [1], the cosmic microwave background 
(CMB) [2] and the large scale structure (LSS) [3] all indicate that the universe is undergoing an accelerating 
expansion. This implies the existence of a mysterious component, named dark energy, which has negative 
pressure and takes the largest proportion of the total density in the present universe. In the last decade, lots 
of efforts have been made to understand dark energy [4], yet we still know litde about its nature. 

In this paper we discuss the possible direct, non-gravitational interaction between the dark sectors in the 
framework of the holographic dark energy model, which is a quantum gravity approach to the dark energy 
problem [5-7]. In this model, the vacuum energy is viewed as dark energy, and is related to the event horizon 
of the universe when we require that the zero-point energy of the system should not exceed the mass of a 
black hole with the same size. In this way, we have the holographic dark energy (HDE) density [8] 



where c is a dimensionless model parameter, which can only be determined by obsei-vations, Mpi is the 
reduced Planck mass, and /?/, is the future event horizon of the universe, defined as 



The HDE model has been proved to be a competitive and promising dark energy candidate. It can theoret- 
ically explain the coincidence problem [8], and is proven to be perturbational stable [9]. It is also favored 
by the observational data [10]. For more studies on the HDE model, see, e.g., [1 1-14]. 

The HDE model with some interaction between dark energy and dark matter (hereafter IHDE model) 
was firstly studied by Wang et al. in [15]. If dark energy interacts with cold dark matter, the continuity 
equations for them are pde + 3H(fide + Pde) = -Q and pdm + ^Hpdm = Q, where Q phenomenologically 
describes the interaction. The interaction between the dark sectors in the HDE model has been extensively 
studied in, e.g., [16]. It was found that the introduction of interaction may not only alleviate the cosmic 
coincidence problem, but also help to avoid the future big-rip singularity. 

There ai^e various choices for the foniis of Q. The most common choice is 



where F is a dimensionless constant, and p is taken to be the density of dark energy, dark matter, or the sum 
of them. These models are mathematically simple, so they are useful for phenomenology. However, it is 
difficult to see how they can emerge from a physical description of dark sector interaction. It is expected 



Pde 



(1) 




(2) 



Q = THp, 



(3) 
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that the interaction is determined by the local properties of the dark sectors, i.e., p^,,, and p^e, but it is hard to 
understand why the interaction term must be proportional to the Hubble expansion rate H. A more natural 
form of the interaction was proposed by Ma et al. [17], where the interaction term takes the form 

Q^PlA- (4) 

In this description the interaction temi only depends on the local energy densities of the dark sectors, and 
is thus more physically plausible. When a = 1 and yS = 0, the interaction term Q oc p^,,, is the exact form 
we expected in the case of the dark matter decay. Similarly, the case of a = and y6 = 1 corresponds to the 
dark energy decay. Moreover, in more complicated forms of interaction, like scattering, one may expect the 
existence of both pd,n and pde- Thus, Eq. (4) seems a more natural and physically plausible form to describe 
the interaction. 

In this paper we revisit the interaction between the HDE and dark matter by considering the interaction 
with the form Q oc p^^^p^^. To perform an overall analysis, we will take into consideration as many factors as 
possible. Different from many works in which the contribution of baryons is ignored, we will consider both 
baryon and dark matter components in the interacting models. Moreover, we will also discuss the models 
in a non-flat universe. As argued in Ref. [18], the studies of dark energy, and in particular, of observational 
data, should include Q^o as a free parameter to be fitted alongside the w(z) parameters. Another reason why 
we take the spatial curvature into consideration is the possible correlation between the curvature and the 
interaction. As pointed out in Ref. [19], compared with the flat universe, a much stronger interaction can be 
allowed by the data in a non-flat universe. 

There are many interesting issues worth investigating in the IHDE models. For example, the properties 
of the interaction term, the direction of the energy flow, the evolution of the dark matter and dark energy 
densities, the fate of the universe, and so on. In the IHDE model with Q oc p^^^p^^, some of these issues have 
been discussed in [17]. In this paper, we will perform a more comprehensive exploration of these issues. 

This paper is organized as follows. In Sec. II, we derive the basic equations for the IHDE models. In 
Sec. Ill, we introduce the methodology and data used in this work. In Sec. IV, we discuss the cosmological 
interpretations of the three IHDE models, including two simple models with Q oc pj„,, Q oc p^g and one 
general model with a and /3 treated as free parameters. We firstly study the cosmological constraints on these 
models, and then discuss the fate of the universe in the models. At last, we give some concluding remarks 
in Sec. V. In this work, we assume today's scale factor ao = 1, so the redshift z satisfies z - 1/a - 1; the 
subscript "0" always indicates the present value of the corresponding quantity, and the unit with c - h = \ 
is used. 
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II. INTERACTING HOLOGRAPHIC DARK ENERGY MODEL IN A NON-FLAT UNIVERSE 

In this section, we give the basic equations for the IHDE model in a non-flat universe. 

A. Friedmann equations in a non-flat universe 

In a spatially non-flat Friedmann-Robertson-Walker universe, the Friedmann equation can be written as 

3Mp;//^ = pdm +Ph+ Pr +Pk+ Pde, (5) 

where pk = -^Mpj-^ is the effective energy density of the curvature component. For convenience, we define 
the fractional energy densities of the various components, i.e., 

- -.T ^ = —, ^de = — , ^dm = = —, = —, (6) 

H^a^ Pc Pc Pc Pc Pc 

where pc - SM'^jH^ is the critical density of the universe. The subscripts, k, de, dm, b and r, represent 
curvature, dark energy, dark matter, baryon and radiation, respectively. By definition, we have 

^de + ^dm + Q.b + ^r + ^k = '^- (7) 

With the existence of interaction between the dark sectors, the energy conservation equations for the 
components in the universe take the forms 

Pdm + 'iHpdm = Q, (8) 
Pde + 3H(pde + Pde) = -Q, (9) 

Pb + 3Hpb - 0, (10) 

Pr + 4Hpr^0, (11) 

Pk + 2Hpt = 0. (12) 

As mentioned above, Q denotes the phenomenological interaction term. Combining Eqs. (8)-(12) together, 
we can obtain the form of pde, 

2 H \ \ 

Pde ^ ~ 3 -Pc- -^Pr + (13) 
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Substituting pde into Eq. (9), we have 



(2- + ^ + 3H)pde + H(pk - pr) - (2- + 3h)p, = -Q. 



'■de 

Dividing the above equation by pc, we get a derivative equation of H and tl^u, 

2{^de - 1)^ + ^de + HOade - 3 + O^t - ^r) = -HQl, 
H 

where we have defined the effective dimensionless quantity for interaction, O/, with the form 

a, = ^ . 

Hiz)Pc 

B. HDE in a non-flat universe 

From the energy density of the HDE, Eq. (1), we have 



Following Ref. [13], in a non-flat universe, the IR cut-off length scale L takes the fonn 

L = ar{t), 

and r{t) satisfies 

r(') dr _ p~ dl_ 
Jo Vl - kr^ Jt a{t)' 

By carrying out the integration, we have 



r{t)^—sm{yk\ — =— sm — -). 



Equation (18) leads to another equation about r{t), namely, 

L c 
r{f) = - 



Combining Eqs. (20) and (21) yields 

c^[k 



Jt a 



arcsm ■ 



Taking derivative of Eq. (22) with respect to t, one can get 



^de ^ ^ ^ ^ _ / ^deH^ k 



2Q.de " H \ c2 a2- 



C. Evolution equations of E(z) and Q^^ (z) in an IHDE scenario 



Combining Eq. (15) with Eq. (23), we eventually obtain the following two equations governing the 
dynamical evolution of the IHDE model in a non-flat universe, 

1 dE{z) a^e 



Eiz) dz l+z 



Q^-Q,-3 + a/ 1 Qrf, 

+ - + a/— 



dz l+z 



^de r\ 1 Q.k- + 
+ ilji + — 



(24) 



(25) 



c2 2 2(1 -Orf,) 

where E{z) = H(z)/Ho is the dimensionless Hubble expansion rate. Equations (24) and (25) can be solved 
numerically and will be used in the data analysis procedure. Notice that we have 

EizY E{zY E{zy 

and the fractional density of dark matter is given by Q-dmiz) - I- ^kiz) - ^deiz) - ^riz) - i^tiz)- The values 
of Ofoo and Qh)> for simplicity, are determined from the 7-yr WMAP observations [23], 

Qto = 0.02253^^ (27) 

Q,-o - Qyoil + O.iniNeff), - 2.469 x I0~^h-^, Ne/j - 3.04, (28) 
where y represents photons, and Ne/f is the efi'ective number of neutrino species. 

D. Models 

In this paper we consider the interaction term with the form Q oc p^^p^^. For concreteness, we express 
the interaction term as 

Q = THo (29) 

which can also be expressed as 2 = THopcoE{zY''"'''^^Q.dmiz)"i^deiz)^, and the parameter F is dimensionless. 

We will investigate three IHDE models in this paper. Firstly, we discuss two simple models with fixed 
a and yS. The first model (hereafter IHDEl) is the case of = 1 and P -0, and thus 

Q = THoPd,n, nj^mdn,iz)/Eiz). (30) 

When F > 0, the energy transfer corresponds to the decay of dark energy into dark matter; vice versa. 
This model is similar to the model proposed in [20], where the same fonii is introduced in the interaction 
between dark matter and quintessence. 



The second model (hereafter IHDE2) considered in this paper is the case of or = and yS = 1 . Cone- 
spondingly, we have 

Q = rHoPde, ni = rQde{z)/E(z). (31) 

Finally, we also consider a more general model (here after IHDE3) where a and fi are treated as free 
parameters. The formula of Q has been given in Eq. (29), while the dimensionless quantity for interaction 
defined previously is 

Qi = rE(zf^"^''^-'ad,„izr^deizf. (32) 



III. METHODOLOGY 



For the IHDE models in a non-flat universe, there are six free parameters: c, Q.dmO, ^ko, r> o; and /?. We 
will constrain them by using the latest observational data. As a comparison, the ACDM model and the HDE 
model with spatial curvature but without interaction (namely, Q^o + but 2 = 0) will also be investigated. 

In this work, we adopt the statistic to estimate the model parameters. For a physical quantity ^ 
with experimentally measured value ^ohs, standard deviation cr^ and theoretically predicted value ^th, ^^X^ 
function takes the form 

= -2 ■ (33) 

The total is the sum of all;t'|s, i.e. 

= 24 ■ (34) 

One can determine the best-fit model parameters by minimizing the total x^- Moreover, by calculating 
Ax^ = ~ A'min' '^^^ determine the 68.3% and the 95.4% CL ranges of a specific model. 

In this work, we determine the best-fit values and the 68.3% and 95.4% CL ranges of the model pa- 
rameters by using the Markov Chain Monte Carlo (MCMC) technique. We modify the publicly available 
CosmoMC package [21] and generate 10^-10^ samples for each set of results presented in this paper. 

For data, we use the Union2.1 SNIa sample [22], the CMB anisotropy data from the 7-yr WMAP ob- 
servations [23], the BAO results from the SDSS DRV [24], 6dFGS [25] and WiggleZ Dark Energy Survey 
[26], and the Hubble constant measurement from the WFC3 on the HST [27]. In the following, we briefly 
describe how these data are included into the^^ analysis. 



7 



A. The SNIa data 



First we start with the SNIa observations. We use the latest Union2.1 sample including 580 SNIa that 
are given in terms of the distance modulus jUobsiZi) [22]. The theoretical distance modulus is defined as 

mizi) = 5 logio DiiZi) + juo, (35) 

where //q = 42.38 - 5 logjg h with h the Hubble constant Hq in units of 100 km/s/Mpc, and the Hubble-free 



luminosity distance Dl = HqcIl is 



Ddz) = -^^M r (36) 

y/Kh^l ^ Jo E{z')' 



where 

sin(;ic), if Q^o < 0, 

Sinn {x) = - x, if Q/to = 0, 

sinh(;c), if fl^^o > 0- 

,1 



The X function for the SNIa data is 



580 



2 _ y' \iiohs{Z.i) - HthiZi)] ^^^^ 



■ 1 cr- 



where jUohsiZi) and cr, are the observed value and the corresponding 68.3% error of distance modulus for each 
supernova, respectively. For convenience, people often analytically marginalize the nuisance parameter fiQ 
(i.e., the Hubble constant Hq) when calculating [28]. 

It should be stressed that Eq. (37) only considers the statistical errors from SNIa, and ignores the sys- 
tematic errors from SNIa. To include the effect of systematic errors into our analysis, we will follow the 
prescription for using the Union2.1 compilation provided in [29]. The key of this prescription is a 580 x 580 
covariance matrix, Cs n, which captures the systematic errors from SNIa (This covariance matrix with sys- 
tematics can be downloaded from [29]). Utilizing Csn, we can calculate the following quantities 

A = itif^ - i//')iCs],)^jifif^ - f,f), (38) 

580 

B = Y,iC-s],)iMf' - Hf), (39) 

i=l 

580 

C - 2 ^^'sN^'j' (40) 

and the function for the SNIa data is [29] 

xIm-A-^. (41) 
Different from Eq. (37), this formula includes the effect of systematic errors from SNIa. 



B. The CMB data 



Here we use the "WMAP distance priors" given by the 7-yr WMAP observations [23]. The distance 
priors include the "acoustic scale" the "shift pai^ameter" R, and the redshift of the decoupling epoch of 
photons z*- The acoustic scale I a, which represents the CMB multipole corresponding to the location of the 
acoustic peak, is defined as [23] 

/a = (1+z*) — — - . (42) 
Here Da(z) is the proper angular- diameter distance, given by 

DA{z) = dL{z)Kl+zf, (43) 



1 r uu 

^siz) = — „ „ z , (44) 



and Vsiz) is the comoving sound horizon size, given by 

^1/(1+-^) da_ 

a^H{a) VI + (3aM)/4f2^o)a ' 
where Q^o and Q.yQ are the present baryon and photon density parameters, respectively. As mentioned above, 

we adopt the best-fit values, Q^o = 0.02253/1"^ and Q.yQ ^ 2.469 x lO"^/!"^ (for Tanb = 2.725 K), given by 
the 7-yr WMAP observations [23]. The fitting function of z* was proposed by Hu and Sugiyama [30]: 

- 1048[1 + 0.00124(Qm/z^)"°"^][1 + giiQmoh^)^'] , (45) 

where 

_ 0.0783(Qfoo/^^)'"-^^^ _ 0.560 

" I + 39.5(Qio/^2)0-763 ' - 1 + 21.1(Qio/z2)i **i ' ^ ^ 

Here the subscript m denotes the matter component, i.e., Q,„ = Cldm + ^b- In addition, the shift parameter R 
is defined as [31] 



Riz,) = ^a^oH^il + z.)Da{z*) . (47) 

This parameter has been widely used to constrain various cosmological models [32]. 
As shown in [23], the^^ function of the CMB data is 

where x; - {Ia,R,z*) is a vector, and {C^\fg)ij is the inverse covariance matrix. The 7-yr WMAP observa- 
tions [23] have provided the maximum likelihood values: Uiz*) = 302.09, Riz*) = 1.725, and z* = 1091.3. 
The inverse covariance matrix was also given in [23], 



(^CMs) 



^ 2.305 29.698 -1.333 ^ 



29.698 6825.27 -113.180 
-1.333 -113.180 3.414 



(49) 
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C. The BAO data 



In this paper we use the BAO data from the SDSS DR7 [24], the 6dFGS [25], and the WiggleZ Dark 
Energy Survey [26]. In the following, we will describe the BAO distance measurements of these projects, 
and introduce how to add them into the statistics. 

One effective distance measure is Dy(z), which can be obtained from the spherical average [33] 



Dviz) 



(l+zfoliz)- ^ 



1/3 

(50) 



H{z) 

where Da{z) is the proper angular diameter distance. 

For the 6dFGS and SDSS DRV data, we use the quantity d- = rsiZd)/Dviz)- The expression of is given 
in Eq.(44), and Zd denotes the redshift of the drag epoch, whose fitting formula is proposed by Eisenstein 
andHu [34], 

1291(Q,„o/l^)0-251 r . 

^'^= l+0.659(a.o/.^)O.S28 [^^^'(^^o^^) ]' (51) 



where 



by = 0.313(Q,„o/t')'''-^'''[l+0.607(a,„o/z')''"'^], (52) 
b2 = 0.238(Qmo/i^)°-^^^ (53) 

For the SDSS DRV data [24], we write for the BAO data as 

XbAO,SDSS ~ ^Pi(^ BA.O,S DS S^ij^P (^'*) 

where 

Api = pdata _ pdata ^ ^ate ^ Q ygQ^^ ^data ^ ^ata ^ q^^q^^^ (55) 

and the inverse covariance matrix takes the form 



(^BAO,SDSS^ 



' 30124 -17227'' 



(56) 



-17227 86977 

The 6dFGS survey [25] gives a measurement of d^^^^^ - 0.336 + 0.015, so we have 

2 _/ ^o.io6 -0.336 \^ 

^BAOfidFGS ~ \ 015 / ' 

The WiggleZ Dark Energy Survey [26] gives three measurements of the A parameter in different red- 
shifts. The A parameter is defined by [33] 

AM.""""''*^. (58) 

z 
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and thex^ of the WiggleZ Dark Energy Survey takes the form 



where 



^BAO y/iggleZ 

- ^PiiCBAO,Wigglez)U^Pj' (^9) 



Api ^ pf'^ - pdata ^ A^f^ = 0.474, = A'^^ = 0.442, ^ - A^f^ - 0.424, (60) 



and the inverse covariance matrix takes the form 



^ 1040.3 -807.5 336.8 ^ 



BAO,Wigglez) " 



(61) 



-807.5 3720.3 -1551.9 
336.8 -1551.9 2914.9 
The final BAO;^^^ is a combination of the SDSS DR7, the 6dFGS and the WiggleZ BAO x^, i e., 

XbAO -^BA0,SDSS ^BAOfidFGS ^BAO, WiggleZ' (^2) 

D. The Hubble constant data 

The precise measurements of Hq will be helpful to break the degeneracy between it and dark energy 
parameters [35]. When combined with the CMB measurement, it can lead to precise measure of the dark 
energy equation of state (EOS), w [36]. Recently, using the WFC3 on the HST, Riess et al. obtained an 
accurate determination of the Hubble constant [27], 

Hq - 73.8 ± 2.4 km/s/Mpc, (63) 

corresponding to a 3.3% uncertainty. So the^^ of the Hubble constant measurement is 

2 //i-0.738f 

E. The total 

Since the SNIa, CMB, BAO and Hq are effectively independent measurements, we can combine them 
by simply adding together the;if^ functions, i.e., 

xltai =xjN+ xImb + xIao + xl ■ (65) 

IV. RESULTS 



In this section, firstly, we show the cosmological constraints on the three IHDE models, and then we 
discuss the cosmological implications of these models. 
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TABLE I: Fitting results of the models 



Model 




c 


c (95.4% range) 


r 


'4-0 


2 

^ min 


ACDM 


235+0012 






(fixed) 


001 +""05 


550.354 


HDE (Q = 0) 


n 9-2^+0.012 


71+0.10 


0.56 <c< 0.94 


(fixed) 


006+" ""^ 

"■"""-0.007 


549.461 


Q = FHopdm 




69+"'° 


0.54 < c < 0.91 


-0 056+""5' 


015+0011 


548.352 


Q = THopde 


0.237^^1^ 


66+"" 


0.50 <c< 0.90 


-0 073+"'"" 

"■"'-'-0.072 


016+""'- 

"■"^"-0.012 


548.390 


Q - ^H,pi/jp:r' 


73Q+0.0I3 
"•^•^^-0.013 


63+""^ 

"■"■^-0.08 


0.47 <c< 0.92 


-0 014+" "'4 

"■"^^-0.237 


017+0.011 

"•"^ -0.014 


548.298 



A. Cosmological constraints 

In this subsection, we discuss the cosmological constraints on the three IHDE models. The parameter 
space of these IHDE models are explored by using the Union2.I+BAO-i-CMB+//o data. We summarize the 
fitting results in Table I. For comparison, the fitting results of the ACDM and HDE models by using the 
same set of data ai^e also shown. 

The values of x^^j^ ^^'^ models are listed in column 7. Clearly, compared with the ACDM and HDE 
models, the IHDE models lead to a evident reduction X^nin- However, taking the number of parameters 
into account, ACDM and HDE models still provide nice fits to the data. In addition, compai^ed with the two 
simple IHDE models, the complicated model with Q oc p^,,jP^^ does not lead to a significant reduction of 
x\^i,^, implying that in the context of the current observations, such a complicated model is not necessary. 

In column 2-6, the best-fit values and the 68.3% CL {Kx^ - 1) uncertainties of the parameters O.dm0, c, 
F and Q-ko are listed. For the HDE and three IHDE models, we found c < 1 at the 95.4% CL {^x'^- ^ A). 
This result is consistent with the results in the previous works [10, 16]. For the HDE model, c < 1 means 
that the universe will end up with a the big rip, while for the IHDE models, due to the interaction between 
the dark sectors, the big rip may be avoided. We will discuss this issue in the following content. 

For these three IHDE models, we find F < at the 68.3% CL, namely, the energy flow from dark matter 
to dark energy is slightly favored by the data. 

In the left panel of Fig. 1, we plot the contours of the 68.3% and 95.4% CL for the three IHDE models in 
the F-c plane. We can see that the upper bounds on the parameter c are similar, while the lower constraints 
are slightly different from each other. In the IHDE3 model, due to the complexity of the model, the allowed 
parameter space is much larger than those of the IHDEl and IHDE2 models. 

Interestingly, it seems that, compared with the HDE model, the IHDE models slightly favor smaller 
values of c. This can be seen in the column 4 of Table I. To see it more clearly, we also plot the probability 
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FIG. 1: Marginalized probability contours at the 68.3% (^x^ = 2.3) and 95.4% (^x^ = 6.18) CLs in the T-c (for 
the three IHDE models) and Q.d„,Q-c (for the HDE model) planes. The contours for the IHDEl (cygan, labeled as I), 
IHDE2 (orange, labeled as II), and IHDE3 (olive, labeled as HI) models are plotted in the left panel, and the contours 
for the HDE model (magnetic) are plotted in the right panel. The green dashed line denotes c = 1 . Clearly, compared 
with the HDE model, the fitting results of c in the three IHDE models are smaller 

contours in the QdmO-c plane for the HDE model in the right panel of Fig. 1. The 95.4% CL of the HDE 
model intersects with the line c = I, while the contours of the IHDE models are all below the line. 

In the left panel of Fig. 2, we plot the contours of the 68.3% and 95.4% CL for the three IHDE models 
in the T-Qi^ plane. The parameter space of the IHDE3 model is much larger than those of the IHDEl and 
IHDE2 models, and we can see clear degeneracy between Q.kQ and F, consistent with the result of [19]. The 
degeneracy amplifies the range of Qto compared with the HDE model without interaction, as is seen in the 
right panel of Fig. 2, where the contours in the Q.dmQ-^kQ plane for the three IHDE and the HDE models are 
all plotted. 

We are also interested in the constraints on the parameters a and /3, i.e., the interaction forms allowed by 
the data. In our analysis, we find the 95.4% CL constraints are 

a = -5.irl',l\ ;S - 10.98!^^:^^ (66) 

So the allowed regions of a and /3 are very wide. This phenomenon is understandable: the intensity of the 
interaction term Q is almost decided by F. If F is small enough, any large values of a and /3 are tolerable. ^ 
The 68.3% and 95.4% CL contours in the a-/3 plane are plotted in Fig. 3, and we can see strong degeneracy 

' In our analysis we generate about 10' samples for the IHDE3 model. It is expected that the area of the contour could be larger 
if one performs a more extensive analysis. 
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0.23 



FIG. 2: Marginalized probability contours at the 68.3% and 95.4% CLs in the F-Qao (left panel) and Qj/mo-Q/to (right 
panel) planes for the IHDEl (cygan, labeled as I), IHDE2 (orange, labeled as II), IHDE3 (olive, labeled as III), and 
the HDE (magnetic) models. The D.di„o-^ko contours of the three IHDE models are similar to each other, so for the 
IHDE2 and IHDE3 models, we only plot the 95.4% CL contours. 
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FIG. 3: Marginalized probability contours at the 68.3% and 95.4% CLs in the a-/3 planes for the IHDE3 model, 
between a and p. 

B. Constraints on the interaction term 

In this subsection we discuss the constraints on the interaction term in the IHDE models. The fitting 
results of T have been listed in the 5-th column of Table I. We have seen that for all the models we have 
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FIG. 4: Reconstructed evolutions of ItHpdm, 3H{1 + Weff)pde and Q (all divided by pco) at the 95.4% CL, for the 
IHDEl (left panel), IHDE2 (middle panel) and IHDE3 (right panel) models, respectively. 

r < at the 68.3% CL. Here we list the 95.4% CL constraints for the three IHDE models, 

-0.154 < r< 0.053 for IHDEl, (67) 
-0.229 < r< 0.059 forlHDE2, (68) 
-0.509 < r< 0.072 forlHDE3. (69) 

For the F > case, i.e., the direction of the energy flow is from dark energy to dark matter, we find a tight 
constraint |F| < 0.07 at the 95.4% CL for all the three IHDE models. For the F < case, the allowed value 
of F is much larger. At the 95.4% CL, for the IHDEl and IHDE2 models we have |F| < 0.2, while for the 
IHDE3 model we have |F| < 0.5. 

To see the result more vividly, in Fig. 4 we reconstruct the evolution of 3Hp4,„, 3H{1 + Wde)Pde, and Q at 
the 95.4% CL for the three IHDE models. As shown in Eqs. (8) and (9), the term 'iHpdm [or, 3//(l +Wde)Pde\ 
characterizes the influence of the cosmic expansion on the energy density of dark matter (or, dark energy), 
while Q describes the effect of the interaction. As shown in Fig. 4, for all the models Q is much smaller 
compared with 3Hpdm and 3H{1 + Wde)Pde, implying that the influence of interaction on pdm and pde is 
much weaker than the cosmic expansion. It is also interesting to find that the evolution of the interaction 
terms have different behaviors in the three IHDE models, due to the different expressions of Q. In the 
IHDEl model with Q oc pd,„, the absolute value of Q may significantly decrease (along with pdm) during 
the evolution, while in the IHDE2 model the absolute value of Q almost maintains no changes during the 
evolution since Q oc p^g. 

We are interested in the influence of the interaction on the dark matter density, and thus we investigate 
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the quantity Q/OHpdm), which describes the ratio between the changes of pj„, caused by the interaction 
and by the cosmic expansion. Here we list the constraints on this quantity at z = at the 95.4% CL, 

- 0.0514 < ^ < 0.0177 for IHDEl, (70) 

-0.2097 < ^ < 0.0625 for IHDE2, (71) 

-0.4779 < ^ < 0.0105 for IHDE3. (72) 

We see that in the IHDEl model, the present effect of Q is much smaller than the other two IHDE models. 
The large value (47.79%) of the lower bound of this ratio in the IHDE3 model indicates that in this model 
the influence of Q on p^/m can be comparable with that of the cosmic expansion. Figure 4 shows that in the 
IHDEl model the present value of Q is much tightly constrained compared with the other two models, and 
in the IHDE3 model Q could take large negative values when z — > 0. 

Also, since we are interested in how much dark matter is produced or annihilated due to the interaction, 
it is helpful to define another quantity, 

Pdma\..-pdma\. 
A = ^^^^ —. (73) 



Pdma^\ 



III! 



Notice that pdm<^^ is the total energy of dark matter in a unit comoving volume. In the absence of interaction, 
PdmCi^ is conserved, and A = 0. Here piiynO^\.^^. is the "initial" value of the total energy of dark matter in a 
comoving volumn. It can be calculated at the early epoch, e.g., z > 5000. Thus, A describes the relative 
change amount of energy in the dark matter component caused by the interaction. 
The results for the three IHDE models are as follows: 

- 15.22% < A < 8.83% for IHDEl, (74) 
-13.3% < A < 5.57% forIHDE2, (75) 
-14.49% < A < 6.57% forIHDE3. (76) 

We find that the results of the three models are not much different from each other. Thus, although the 
evolution of the interaction term is evidently different, the constraint on the overall effect of the interaction 
does not change a lot. Roughly, at the 95.4% CL, the increment amount of the energy of dark matter is 
constrained to be less than 9%, while the decrement amount is constrained to be less than 15%. 



C. The fate of the universe 



In this subsection we discuss the fate of the universe in the IHDE models. The most important reason for 
the introduction of the interaction between the dark sectors in the HDE model is to avoid the future big-rip 
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singularity. For the IHDE model with Q oc p^„jP^,„, this issue has been briefly discussed in [17], where it 
is shown that back reaction efi'ect from quantum corrections cannot prohibit the occurrence of the big rip. 
Here, we will give a more detailed investigation of the fate of the universe in the IHDE models considered 
above. 

Firstly, we will derive some useful formulas as a preparation of our discussion. Then, we will analyze 
the dynamical equations to see whether it is possible to avoid the big rip when c < 1 . Finally, based on the 
constraints of the data, we will numerically solve the equations and reconstruct the evolution of the dark 
sectors to z ^ -1. 

1. The effective EOS of dark sectors 

We can define the effective pressure and eff'ective EOS of dark energy, i.e., 

Peff,de = Pde + Weff,de = Wde + tS~' ^^'^ 

satisfying 

Pde + 3//(pde + Peff) = 0, pde + 3^/(1 + Weff4e)Pde = 0. (78) 

From Eq. (13), using = -{I + z)Hj-^ and ^ = -(1 + z)^^ = -(1 + z)^, we derive the equation 
of state for dark energy, 

2{l+z)dlnE 1 ^ ^^^^ 

3Q.de dz Qde 3 Qde 3 Qde 

Combined with Eq. (24), it follows that 



1 2 Qde „ Qi 

which reduces to the familiar formula Wde = ~\ ~ f "^T^ when 0.1^ = and Q.j = 0. Thus, the effective EOS 
takes the form 

Weff,de = Wde + ^ ,P (81) 



1_2 Ode 

■3 3 V c2 



= ----\l^ + ^k- (82) 



This result is very interesting: Weff,de does not explicitly contain the interaction term. However, we should 
keep in mind that the evolutions of Qdeiz) and E{z) are determined by the differential equations (24) and 
(25) where the interaction Q (or, Q/) is involved. 
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It is much easier to derive the effective EOS of dark matter. From Eq. (9), it follows directly that 

Q 

Weff,dm = = (83) 

- -lr£(z)2(-/^)-3Q-if/^^. (84) 

To investigate the fate of the universe, we are also interested in the effective EOS of the total components 
in the universe, i.e.. 



Weffjot 



WeffJ^i, (85) 



where / represents de, dm, b, r and k. In the future epoch, only the dark energy and dark matter components 
are important. 



2. Analysis of the dynamical equations in the z — > — 1 region 

Now we make some general investigations to see whether the big rip can be avoided when c < 1 in the 
IHDE models. Equation (82) shows that the effective EOS of HDE is only determined by Q^g and c (we 
neglect 0.]^ which is less important in the future). In the case of c < 1, once — > 1 there will always 
be Weff4e < -1 and Weff\tot < -1> and the fully dominated dark energy will drive the universe to a big 
rip. Thus, to avoid the big rip, we must evade n^/g — > 1 in the future. Also, to avoid p^e from infinitely 
increasing, we must require F > 0, i.e., the direction of energy flow is from dark energy to dark matter. 

Here we list the simplified equations of motion in the region z ^ -1. Neglecting the radiation, baryon 
and curvature components, we have 

+ , (86) 



E{z) dz l+z\ l^de 



dQ.de 2Q^/e(l - fide) i ^Qde }_ _ ^/ i ^g^-j 



dz l+z \ c 2 2(1 - Qde) 

Notice that unlike the cases of 2 oc Hpd,,, or Q oc Hpde, the two equations are coupled, and thus the situation 
is much more complicated. However, we find that these equations have a stable point 

Q.de = c^, Qi = 3(1 - c^), (88) 



which leads to 



1 2 yjQde Q-I 
^eff4e = -3 - = -1, W,ffM.n = = (89) 
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FIG. 5: Evolutions of E, Q,/^, Qf/,„, w^jfj^ and Weff\dm along with the scale factor a. Here we take typical values 
of QdmO - 0.24, Ojo = 0, /i = 0.73, and others are denoted in the panel. 

and 

^=0, ^=0. (90) 
dz dz 

This is exactly a de Sitter solution. See also [37] for a relevant study about the avoidance of big rip in the 
holographic dark energy model, where the extra dimension effect is considered and the final state of the 
universe is also a de Sitter solution. 

We find that this situation can really happen when y6 is large, and the interaction is strong enough to pull 
down Qf/e from 1. As an example, in Fig. 5 we demonstrate how the quantities evolve along with the scale 
factor a in the case of /? = 10. A sophisticated discussion on the evolution behaviors of the quantities in the 
whole range of a would be tedious and unnecessary. Thus, instead, we just focus on the epoch of a > 10 
(the region to the right of the thick dashed line). This would be enough for us to understand how it happens. 
In this epoch, firstly, E{z) increases due to the dominated dark energy component with Weff\de < -1, while 
the strong interaction pulls O^g down from 1. So along with a, Q^e decreases and correspondingly Qj,,, 
increases. As a result, Wgff^de increases from values less than -1, while the effective EOS of dark matter, 
Weff4m - 3n^' which is highly negative due to a small Q.dm, also goes up. When both Wde,eff and Wdm,eff 
approach -1, the system reaches the stable point, and E{z) stops increasing. The final state of the universe 
is de Sitter-like. 

Note that an essential condition that such a procedure can happen is that the interaction is strong enough 
to pull Qde down from 1 . Thus, a large /3 is expected. To see the dependence on the values of j3, in Fig. 6 
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FIG. 6: Evolutions of Udm (left panel), 0.^^ (left panel), We/f.de (middle panel) and Weffjm (right panel) along with the 
scale factor a for three sets of a and p. The same as Fig. 5, we take values of Q^/,„o — 0.24, Q.k(i - and li - 0.73. 

we plot the evolution of Q^/g, Q.dm, Weff,de and Weff^dm in large a regions for different values of a and /3. For 
the case of a = and P - \ (the blue line), the interaction is not strong enough to avoid Q.de — > 1, so the big 
rip still happens. For the case of a = and fi = 3 (the green dotted line), the interaction is strong enough, 
and so we see that Q.de is pulled down from 1, and the system evolves to a stable point when a > 10000. 
For the case of a = -1 and /3 = 7 (the red dashed line), the interaction is so strong that the system quickly 
approaches the stable point at a < 1000. 

In the following, we will numerically solve the equations for the IHDE models in a large a region in 
the parameter space constrained by the data, and see whether the de Sitter solution can be achieved in these 
models. 



3. The IHDEl and IHDE2 models 



In the top panels of Fig. 7 we show the reconstructed evolution of Qde along with z for the IHDEl and 
IHDE2 models. Clearly, at the 95.4% CL we see Qjg 1 when z ^ -1, and we have Wef/.de < -1 since 
c < 1. So these models will not help us to avoid the big rip. 

This result is understandable. The fate of the universe in the IHDE model with the interaction term 

Q - 3bHpde (91) 
was discussed in [9] and it was shown that the condition to prevent the big rip is 

b>c-^ - 1. (92) 
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FIG. 7: The evolutions of Q^^ (top panels) and Weff,de (bottom panels) along with the redshift z at the 68.3% and 
95.4% CLs, for the IHDEl (left panels) and IHDE2 (right panels) models. The 95.4% CL of HDE model is plotted in 
dashed lines for a comparison. 

In the IHDEl and IHDE2 models, analogously, we define the "effective" coupling, 

beff = — (93) 

and we have 

^ _ ^^dm ^ _ r (94") 

for the two models, respectively. Note that once E{z) becomes increasing in the future, be// will be sup- 
pressed to zero, so the condition (92) cannot be satisfied. 

The reconstructed Weff\de(z) for the two models are shown in the lower panels of Fig. 7. As a comparison, 
the 95.4% CL evolution of the HDE model is also plotted (the black dashed line). We see that for these 
models Wg/f^deiz -1) < -1, and thus the big rip will happen. Also, the values of Weffje in the IHDE 
models are more negative than that of the HDE model, since in the two IHDE models c is constrained to 
be smaller values. Thus, instead of helping us to avoid the big rip, the situation is even worse in these two 
models. 
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FIG. 8; The evolutions of Side (gray, left panel), D.dm (green, left panel), Weff,de (gray, right panel) and Weffjot (orange, 
right panel) along with the scale factor a for the IHDE3 model. The 68.3% and 95.4% CLs are plotted. 

4. The IHDE3 model 

We have seen in Fig. 6 that a large yS can drive the system to a stable point and avoid the big rip. As 
shown in Eq. (66) and Fig. 3, in our fitting results we found that large values of /? are allowed, so we expect 
that the big rip can be avoided in the IHDE3 model. 

In Fig. 8 we plot the evolutions of Q-de (gray, left panel), Orf„, (green, left panel), Weff,de (gray, right 
panel) and Weff,tot (orange, right panel) along with the scale factor a for the IHDE3 model. From the left 
panel, we see that, at the 68.3% CL Q.de — > 1 is present in the future due to a negative F (see Table I), but 
the interaction can prevent Q^e from approaching 1 at the lai^ge a region at the 95.4% CL. The right panel 
shows that the stable point (corresponds to Wef/^de - Wgffjot - -1) can be accomplished when a > 1000 at 
the 95.4% CL. Thus, the IHDE3 model may help to avoid the big rip singularity. 

It would be worthy further investigating what kind of samples in our MC analysis may achieve the 
de Sitter solution. Of 10^ samples generated for the IHDE3 model there are about 4.7 million samples 
satisfying Ax^ < 4. Within them, we find about 18,000 samples satisfying F > 0, and finally 709 samples 
having de Sitter solution in the far future. In Fig. 9, these samples are plotted in the a-yS (top panels) and 
F-yS (bottom panels) planes. Clearly, we see that to achieve a de Sitter solution, a large yS is required (all 
these samples have yS > 1.5). The bottom right panel of the figure shows that, at the 68.3% {Ax^ < 1) CL, 
due to a negative F, the big rip cannot be avoided. 
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FIG. 9: b^x^ - 1 (olive dotted) and Ay^ = 4 (green dashed) contours in the a-fi (top panels) and F-jS (bottom panels) 
planes. Orange dots are the 709 samples which can achieve a de Sitter solution in the future and thus avoid the big rip. 
In the left panels the whole contours are plotted, while the right panels, as the zoom in of the dashed orange region of 
the left panels, show these samples in detail, ji - 1 .5 is plotted in black dotted lines. 

V. CONCLUDING REMARKS 

We have investigated the IHDE models with Q oc P^,,,/^^- Three IHDE models, including the IHDEl 
model with Q = THopd,„, the IHDE2 model with Q = YHopde, and the IHDE3 model with a and p running 
freely, were investigated. 

By using the Union2.1-i-BAO+CMB-i-//o data, we placed cosmological constraints on these models. We 
found that a negative F, i.e., energy flow from dark matter to dark energy, is slightly favored by the data, 
although the HDE model with F = still lies in the 95.4% CL region. For all the models, we get c < 1 at 
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the 95.4% CL. The values of c in the IHDE models are smaller than that in the HDE model. 

We showed that the interaction has different properties in the thi^ee models. Compared with the cosmic 
expansion, the effect of interaction on the evolution of p^m and pde is smaller. We also put a constraint on the 
total amount of the energy change in the dark matter component. At the 95.4% CL, the increment amount 
of the dark matter is constrained to be less than 9%, while the decrement amount is constrained to be less 
than 15%. We found that the constraint basically does not depend on the forms of the interaction term. 

Furthermore, we discussed the fate of the universe by investigating the dynamical equations of the IHDE 
models. We find that the equations may give a de Sitter solution at z -1, with the effective equation of 
state for dark energy and dark matter being -1. When confronted with data, we show that this solution 
cannot be accomplished in the IHDEl and IHDE2 models. Rather, in these two models since the data favor 
a smaller c, the big rip is even more severe than the HDE model. In the IHDE3 model, we show that such a 
solution can be achieved for a large yS, and the big rip may be avoided at the 95.4% CL. 

Finally, there are still some issues not covered in our paper, i.e., the precise condition for obtaining the 
de Sitter solution, the coincidence problem, and so on. These issues all deserve further investigations. 
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